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This work develops a new generalized technique for determining the static and dynamic properties 
of any non-collinear magnetic system. By rotating the spin operators into the local spin reference 
frame, we evaluate the zeroth, first, and second order terms in a Holstein-Primakoff expansion, and 
through a Green's functions approach, we determine the structure factor intensities for the spin-wave 
frequencies. To demonstrate this technique, we examine the spin-wave dynamics of the generalized 
Villain model with a varying interchain interaction. The new interchain coupling expands the overall 
phase diagram with the realization of two non-equivalent canted spin configurations. The rotational 
Holstein-Primakoff expansion provides both analytical and numerical results for the spin dynamics 
and intensities of these phases. 

I. INTRODUCTION 

Disordered and frustrated magnetic systems have provided condensed-matter physics with a range of complicated 
problems [l|. Magnetic frustration is typically induced by competing nearest and next-nearest neighbor antiferro- 
magnctic (AF) interactions JlJJd • To minimize the energy, the local moments frequently rotate into a non-collincar 
ground-state configuration [2|, |3( . Theoretical and experimental efforts have worked to identify and characterize non- 
collincar magnetic systems. Quantum Monte Carlo and Hartree-Fock calculations have been used to describe canted 
antiferromagnetic (CAF) phases in double quantum dots and bi-layer quantum Hall systems H [B[, while ab initio 
methods have been used to examine non-collinearity in magnetic atomic chains . Through neutron-diffraction stud- 
ies, layered borocarbide systems -RB2C (R = Dy, Ho, and Er) have demonstrated conventional and unconventional 
magnetic correlations Q- Other systems like the cuprate, ruthenate, and manganite systems (where the competition 
between AF and ferromagnetic (FM) order has been examined in great depth [T3, [Til, flo| ) also display non-collinear 
characteristics @, i, OH E, E Hi 111 El- To help understand non-collinear magnetic systems, we present a new 
technique for modeling the static and dynamic properties of canted local moments in any periodic system. 

After rotating the spin operators into the local reference frame at each site, we use a Holstein-Primakoff (HP) 
expansion to determine the classical spin energy and the spin-wave (SW) frequencies. By solving the cquations-of- 
motion for coupled Green's functions, we determine the structure factor (SF) intensities for any eigenfrequency of the 
system. While a rotation into the local reference frame of the moments has been employed before to describe the SW 
modes of spin glasses and other non-collinear two-dimensional AFs [H, [l?], [H, HH , the general technique introduced 
in this paper can easily describe any magnetic system, with arbitrary exchange coupling and single-ion anisotropy. 
The inputs into this method are the interaction parameters and moment angles; the outputs are the SW frequencies 
and intensities. 

To demonstrate this general technique, we investigate the generalized Villain model with an added interchain 
coupling. In recent years, the generalized Villain model (GVM) has been frequently used to test approaches to 
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FIG. 1: (a) The generalized Villain model (GVM') with Heisenberg interactions J, 7J, and —r/J. (b) The GVM' with anti- 
parallel spin configuration (SCI), where the local moments of sites a and b reside in the xz-plane and are described by the 
angles (8 a ,6b)- (c) The GVM' with parallel spin configuration (SC2), where the local moments of sites a and b reside in the 
zz-plane and are canted by the angle a. Note that SC2 with a = 7r/2 is rotationally equivalent to SCI with (0, 7r). 
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study magnetically frustrated systems and to help understand more complicated magnetic systems similar to those 
mentioned above [TH, HJ, • In the GVM ; chains of FM interactions J and AF interactions — rjj along the x-axis 
are coupled together in the y direction by J [HI, I20I |2lL [2^ |. as illustrated in Fig. [lja). The original Villain model 
assumed that rj = 1, which denotes full frustration This model was generalized to rj 7^ 1 by Berge et al. 0], 
accounting for coupling of different magnitudes within the FM and AF chains. Within the GVM, a CAF phase is 
stabilized for rj > r\ c . The critical value increases with applied magnetic field from a value of i)c=\ hi zero field. In 
1992, Saslow and Erwin [l8[ employed a HP expansion to numerically examine the SW and SF intensities for this 
generalized model. 

In this paper, the GVM is extended further by introducing a variable y-axis coupling 7 J, as shown in Fig. [ Ha) . We 
call this model the GVM'. Therefore, the GVM' with 7 = 1 reduces to the GVM discussed previously pi MM [2^ . 
As in the GVM, all interactions are confined to the xy plane but the magnetic is applied along the z axis. Through 
an examination of the classical energies, we demonstrate that there exists two separate canted spin configurations 
throughout {77,7} phase space. Figures QJb) and (c) show the possible canted spin configurations within the GVM'. 
In Fig. Hlb), the local moments are canted in the zz-plane by angles 9 a and 9^. In this co nfig uration, the spins 
projected onto the xy plane are antiparallel on sites a\ and a-i as well as on sites b\ and &2 jl8l. I20I l2ll . H2]. This 
canted spin configuration occurs in the GVM with 7 = 1. With the introduction of the variable 7, a new spin 
configuration arises: the corresponding projected moments are parallel with moment angles (a), as shown in Fig. 
[He). For notational convenience, we call these SCI (Spin Configuration 1) for the anti-parallel case (Fig. [Ijb)) and 
SC2 (Spin Configuration 2) for the parallel case (Fig. HJc)). 

Through an examination of the classical limit, we determine the phase diagram within the parameter space {77,7}. 
The GVM' is found to support three phases (FM, SCI, and SC2). Using the rotational Holstein-Primakoff expansion, 
the SW frequencies are determined analytically for all phases. The SF intensities for the FM and zero-field SC2 phases 
are determined exactly; the rest are solved numerically. 



II. THE GENERAL ROTATION MODEL 

As described in Ref. [2l|. the Hamiltonian for canted magnetic systems can be simplified by rotating into the 
reference frame for each moment: S,; = C/jSi, where ]J_ i is the unitary rotation matrix for site i (discussed in Appendix 
A). In the classical limit, Si points along its local z-axis. The general Hamiltonian is given by 



h = - \ E Jifr- Sj - E - 5 E s * 
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where Si are the local moments for site i , is the interaction between sites i and j, Di is the single-ion anisotropy, 
and B is the applied magnetic field. It should be noted that only depends on the moment angles of the spins at 
site i. The canting of a local moment can be described as a rotation by 9 in the xz-plane, with another rotation by %j) 
in the ccy-plane. Therefore, each local moment can be described by the Euler angles 9 and ip (23|. The Hamiltonian is 
expanded in powers of 1/\J~S about the classical or high-spin limit: H=Eo+Hi+H2 + ---- Within the HP formalism, 
the spin operators in the local reference frames become: Si Z — S — a\cn, Si+ = \/2Sdi and Si- = \/2Sa\ Q. The 
zcroth-ordcr Eq term corresponds to the classical energy and the second-order term H2 describes the dynamics of 
non-interacting SWs. The first-order term Hi vanishes when the local moments minimize the classical energy Eo for 
a specific interaction pair rj and 7. Each term, up to second order, is discussed futher below. Higher order terms 
correspond to SW interactions that are unimportant at low temperatures and for small 1/S. 

A. Zeroth Order: Classical Energy 

From the above Hamiltonian, the zeroth-order terms describe the classical energy and can be written as 

E » = ~\ E 'ArW-li - E D i Sfcos(9. l ) 2 --B'E Sfcos(^), (2) 

i.j i i 
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where Si is the moment magnitude for site i, B' = B/ JS, and F\\ is a rotation coefficient given by the angle rotation 
matrix in Appendix A. It only depends on the angles of the spins on sites i and j. From an examination of the global 
minima, the appropriate spin configuration can be determined by minimizing this classical energy. 



B. First Order: Linear Terms 



The HP expansion produces terms that are linear with respect to the creation and annihilation operators. Therefore, 
the first-order Hamiltonian is given by 

#i = - E Jii {^-^( F i a \+ F i* a i) + —7g= i F 2 a ]+ F 2* a j)) +E \/f- {B'+2D i sm{6 i ) cos(fc)) (oj+o,), (3) 

where F[ 3 = + iF l y { and = Fg + iF%. Here, F l x {, F%, F l y {, and F« are rotation coefficients that depend 
only on the moment angles for the interacting spins on sites i and j (described in Appendix A). The linear terms 
correspond to the creation and annihilation of SW's from the vacuum. Assuming the system is in a proper ground 
state with the angles that minimize the energy, the nonphysical first-order coefficients of a; and a\ must vanish for 
each spin site. 



C. Second Order: Spin Dynamics and Structure Factor Intensities 



The second-order terms in the HP expansion describe the spin dynamics. We Fourier transform the spin operators 
by: a k = l/Nj^'f'' e _jkRi ai and a k r ^ = l/Nj^i^ e lkRi a|, where the sums are restricted to sub-lattice (SL) r. For 
each SL of the system, the moment angles are the same. The second-order terms of the generalized Hamiltonian can 
then be written as 

r,s u,k 

+4 S)t 4 S) )} - \ E D r S r (sH8 r f(a^a^ + + a«a™ + ~ ^os^fa^) (4 ) 

r.k 



-#£&a>s(0 r )4 r)t 4 r) > 

(u) . . . (u) 

where k is the momentum vector. Here, z)- s is the number of SL s sites coupled by the interaction Ji- S to a site on 
SL r, and u denotes the multiple possible interactions from SL r to SL s. For example, a FM with a single SL can 
have both nearest-neighbor and next nearest-neighbor interactions with u = 1 and 2. We have also defined r^ 5 = 
1/zr" J2d e~ ik ' d< } with d' u ) = ~Rj — Ri where R; on SL r and Rj on SL s are coupled by the exchange Jr" . Note 
that r^g' 1 = 1 and r^**""' 1 = r[^ u )* = p* r ( u )* p ma Hy 5 Qrs anc j qts are r0 ^ a ti n coefficients that depend only on the 
moment angles for the specific SL (described in Appendix A). 

To determine the SW frequencies oj kl we solve the equation-of- motion for the vectors Vk = 
[a^\ a^^], which may be written in terms of the 2s f x 2s t matrix M(k) as 

idv k /dt = - [H 2 ,v k ] = M(k)v k , (5) 

where s t is the number of spin SLs. The SW frequencies are then determined from the condition Dct[A/(k) — lo^I] = 
0, where only positive frequencies are retained. 

The structure factor in a magnetic system describes the intensity expected from experiment for SW modes with 
resolution- limited width [24j. In the case of a standard FM with identical nearest-neighbor interactions, the structure 
factor is constant throughout fc-space. However, as the spins cant, the wave-vector dependence becomes important. 

Local stability in a magnetic system requires two conditions: I) All SW frequencies must be real for every k and 2) 

the SW weights must be positive. The weights are given by the coefficients of the delta functions in the spin-spin 
correlation function 



S(k,cj) = i[s+'-(k,w) + S-'+(k,a;)] + S z ' z (k,co) = ^ W^S(u - u^), (6) 
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where u) is the eigenfrequency for the spin-wave [24j, |25[ and the sum over t means a sum over all SW modes. The 
total number of transverse or longitudinal SW modes in the first Brillouin zone equals the number of magnetic SLs 

St- 

Generally, the spin-spin correlation function is 

S°"^k,u) = 1 fdt e--*^e*-( R J- R ')(Sf(0)Sf(t)>, (7) 

i,3 

where a, /? = +, — , and z 0, HE] The transverse terms (a, (3 = +, — ) correspond to .S +,_ (k, uj) and 5 _,+ (k, uj) 
in 5(k, w), while the longitudinal term S z,z (k, uj) with a = (3 = z is only nonzero when the system is canted. 

Here, (Sf(0)S?(*)) = ( (U^S^O)) (Uj 1 S j (t)) f3 ) is rotated from the local-moment frame to the global frame. By 

expanding and solving for the spin Green's functions, we can write S a ^(k,oj) as 

S a ' (k,u) = -- limlm (Tr (G(k,w + i8)C a!3 )) , (8) 

where 

G(k, u + iS) = f 1 dTe^' T g k (r) , (9) 

u>i = 2lirT, and C™^ 1 are the rotational coefficients for the canted moments, which depend only on the angles 9 and 
ip (described in Appendix A). Here n = 2r -1, 2r for SL r and 

9nm( T ) = -(T T v„ k (r)v mk ), (10) 

where w„k(T") = e H ' 2T v n yeT H ' 2T . Using the cquations-of-motion of v„k(r), the Green's function matrix can be solved 

as 



G nm (k,oj + iS) = 



1 -N 



(uj + iS)l- M(kY 



(11) 



where / is an 2st x 2st identity matrix and N_ is an 2st x 2st matrix defining the commutation relations such that 
K,k, v m ,k] = N nm . 

From the commutation relation, [S^~ , S~] = 2S z Sij, one can determine the net magnetic moment from the sum rule 



/oo 
duS M Qt,w), (12) 
-oo 



where 5 2 (k, u) = i(5 +_ (k, uj)-S~ + (k, lj)). Note that S z (k,Lj) is not the same as the longitudinal spin-spin correlation 
function, S z ' z (k, uj). Since (S z ) does not depend on k, neither does the right-hand side of Eq. (fT2|) . 



III. THE GENERALIZED VILLAIN MODEL: CLASSICAL ENERGY AND PHASE BOUNDARIES 

To demonstrate the rotational technique described above, we examine the GVM' described in Fig. [TJ This model 
neglects anisotropy, but includes an applied magnetic field. The introduction of 7 into the GVM' expands the 
overall interaction phase space, which can illuminate some of the interesting phenomena seen within the cuprate and 
manganitc systems. An overall analysis of the numerical and analytical results will provide detailed information about 
the GVM', as well as a unique look at the nature of the spin phases as one moves through the interaction phase space. 

While the general rotational Hamiltonian describes interactions between spin sites with moment rotations in both 
the xz (9) and xy (ip) planes, the GVM' constrains the moment angles to rotations in the xz-plane, which greatly 
simplifies the technique. Therefore, spin-spin interactions can be denoted as angle pairs (0 a ,#&). This constraint also 
allows the SW frequencies for the GVM' to be determined analytically. 

Through an examination of the classical energy we determine the energy boundaries for each spin configuration. 
From Eq. ([2]), the classical energy of SCI is 



p(SCi) _ JS 



(?7cos(20 6 ) - cos(20 a ) - 2 7 cos(0 a - 9 b ) - B'(cos{6 a ) + cos(0 b )) , (13) 
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FIG. 2: (a) The phase diagram for the GVM' at B = in the interaction space of 7 and r\. The green/black dotted line indicates 
the phase boundary between SCI and SC2 and the blue dashed line shows the sub-space of the GVM' (7 = 1). (b)-(d) Pictorial 
representations of the (b) (7r/2,7r/2) and (c) (-7r/2,7r/2) limit regions (black dash-dotted line) of the SCI phase, and (d) the 
SC2 phase at zero field. 



(SCI) 

Minimizing Eq with respect to 9 a and 9b yields the relations 

sin(20 a ) + 7sin(0 a - b ) + is'sin(0 a ) = 0, 

(14) 

t]sw(29 b ) + 7 sin(0 o - 9 b ) - ^B'sm(6 b ) = 0. 

Assuming that (9 a ,9b) satisfy these criteria, the first-order terms from the Hamiltonian ii{ will vanish. At zero 
field and in the limit of large 77, the local moments cant toward the angles (tan _1 ( 7 /-\/4 — 7 2 ),7r/2). In the limit 
7 — >1, Eqs. (fTS]) and (fT4|) agree with the results of Ref. [2l|, where 9 a and 9b approach angles smaller than 7r/6 and 
7r/2, respectively. With increasing 77 and | 7 |, the zero-field SCI phase approaches one of the planar phases shown in 
Fig. OJb) and (c). 

By linearizing Eq. (fT4|) , we obtain the phase boundary between the SCI and FM phases : 

B' = 2(77-7- 1± vV + ^+l) 2 )- (15) 

In the limit 7 — >1, this gives the relation obtained by Gabay et al. [26|. Solving Eq. (fT5|) for 7, we find the phase 
boundaries 

(4 + B')(4r,-B') 

7 = ± 4(2-2^ + 5') ■ (16) 
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FIG. 3: Phase Diagram for the GVM' at B' = B/JS — 5 in the interaction space of 7 and r\. The green/black dotted line 
indicates the phase boundary between SCI and SC2. The blue dashed line shows the sub-space of the GVM' where 7 = 1. 



Figure [DJa) shows the phase diagram for the zero-field GVM'. Here, the blue dotted line goes along 7 = 1, where the 
critical value r\ c = 1/3, is consistent with previous work Q. 

In zero field, SCI with angles (0,7r) is equivalent to SC2 with a = tt/2. However, due to the absences of anisotropy 
in the GVM', a magnetic field in the z direction immediately produces a spin flop into the xy plane with the spins 
canted towards the z-axis. This configuration corresponds to SC2 with a < tt/2. Figure^d) shows the SC2 phase at 
zero field, where a = tt/2 throughout the region. Equation (fT6|) gives the field dependence of the SC1/FM boundary. 

To find the field-dependent boundaries for the SC2 phase, we examine the classical energy 



E, 



(SC2) 



JS 2 



rj - 1 - 27Cos(2a) - 2B'cos(a) 



(17) 



At zero field, SC2 has the same energy as SCI. However, SC2 has the lower energy when a magnetic field is applied, 
so it does not follow the same field dependence as SCI. By minimizing the energy, we obtain the SC2 angle 



a = tan 



B' 



4- 



(18) 



which depends only on B' and 7. 
phase boundary is given by 



Consequently, the SC2/FM phase boundary is given by B' = —4j and the SC1/SC2 



7 



2'/ 



(r?-l) 



(19) 



which is shown by the green/black dashed lines in Figs. [2Ja) and [31 This reveals that the SC1/SC2 boundary does 
not depend on field. As the magnetic field increases, the FM phase expands along this boundary, creating a triple 
point at 77 = B' j(B' + 8) and 7 = —B'/4. The expansion of the FM phase space with field is clearly illustrated in 
Fig. H where B' = 5. 



IV. GENERALIZED VILLAIN MODEL: SPIN- WAVES AND STRUCTURE FACTORS 



Using the methods described in section 2, we apply this new technique to the spin dynamics of the GVM' in the 
three regions of the phase diagram (FM, SCI, and SC2). Since the environments surrounding the moments at sites 
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a\ and a 2 are equivalent, they can be considered part of the same sublattice. Similarly, for sites b\ and b 2 - The SW 
frequencies can be determined analytically for all phases. The SF intensities will be obtained for a more complete 
picture of the system. In the FM and zero-field SC2 phases, we solve the SF intensities analytically. However, due to 
the added complexity of the nonzero-field SC2 and SCI phases, those intensities must be determined numerically. 
The second order GVM' Hamiltonian can be written as 

h 2 = jsj: (4 r)t 4 s) 4 r ' s) + + &* )f rt r ' s) ) , (20) 

k,r,s 

where and B^ are coefficients that describe the interactions while r, s = a or b represent the two SLs 

[2l| . The resulting SW frequencies can be expressed analytically using the modified coefficients for the different spin 
configurations, with expressions given in Appendix B. 



A. The FM Phase 



Both the SW frequencies and SF intensities of the FM phase can be obtained analytically. The SW frequencies are 
described by 

^ k = B + JS(2 7 + (77 - 1) (cos (k x ) - 1) + E±) , (21) 

where 

Rf u = ±^J (77 + 1) 2 (1 - cos(fc x )) 2 + 4 7 2 cos(fc y ) 2 (22) 

and the lattice constant a has been set to 1. In the limit k x = 0, u>u. = B + 2775(1 ± cos (k y ) ) only depends on 7 J, 
which is consistent with previous results [H, Hi], H2| ■ As a function of 7, 77, and k, the SF intensities are given by 

Wk = Rf k + 2 7 cos(k y ) _ (23) 

This describes the SF intensities for both the low (+) and high (-) frequency SW. 

Figure [4] shows S(k x ,rj) with 7 = 3 and k y = for the high (a) and low (b) SW modes of the FM phase. If 7 
is held constant, the plot forms a saddle, where the intensity at 77 = -1 is the maximum throughout k x /ir. If 7 > 
1, then this maximum is constant over k Xl which is consistent with the standard two-dimensional FM. When 7 < 
1, the maximum displays non- linear behavior as k x goes to 7r. The sharpness of the saddle depends on 7. As 7 
approaches 0, the saddle sharpens to a delta function at 77 = -1, while increasing 7 flattens the saddle. The mode 
with the maximum intensity is also determined by 7. If 7 > 0, then the low-frequency mode dominates the intensity. 
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FIG. 5: (a) Transverse low- and high-frequencies SW modes at B' = and 77 = 4 for 7 = 1...4 as function of (k x /n,k y /it). 
Their corresponding SF intensities over the same values are shown below. The higher-frequency modes have a lower intensity 
in the k x direction, (b) Longitudinal modes and intensities activated by the canting of the local moment. 



The opposite is true for negative 7. As rj departs from -1 the SF intensity becomes more dependent on k x /-K and 
the intensity decreases. Therefore, even though the system is a two-dimensional FM, the non-equivalent interactions 
create a unique intensity pattern. When the intensities are summed over the high and low frequencies, the total 
intensity is constant throughout k. This demonstrates that the complex interactions distribute the intensity between 
the two SW branches, but the total intensity still remains constant as in the standard ferromagnet. 



B. The SCI Phase 



The SCI region is described by angles (9 a ,6b) that are given by Eq. (fl~4")) . The SW frequency coefficients for the 
GVM are given by Eq. (|B1[) in Appendix B. In this section, we present the SW frequencies and SF intensities for 
varying 7, 77, and B' to give a representative cross-section of the CAF region. Due to the complexity of the SF 
intensities, they have been determined numerically. 

The region 7=1 has been examined in great detail in Refs. [l8[ and [2l|, where the SW intensities increase 
dramatically as [k x / V 7 k y / V) approaches (1,0). A more moderate shift in the intensity is seen along (0,k y /n). As 
shown below, similar features are observed as 7 is increased. 

Fig. [5] shows the low and high SW frequencies and SF intensities for the longitudinal and transverse modes at zero 
field with ?y = 4 and 7 = 1...4. The lower SW frequencies tend to decrease in the (k x /Tr,0) direction and increase in 
the (0,ky/Tr) direction. This difference aries because the exchange 7J is along the y direction. As 7 increases, the 
shift in the intensities along the (Q,k y /-K) direction become more pronounced. As 7 decreases, the SW intensity along 
the (0,k y /ir) direction disappears as the FM and AF chains become decoupled. 

As shown in Fig. [U^a), the corners of the zero-field phase diagram consist of two planar regions sketched in Fig. [2] 
(b) and (c). A magnetic field causes the spins in both regions to cant. In Fig. [51 we investigate the transition from the 
SCI phase into this planar regime. While the SW frequencies do not show much overall difference, the SF intensities 
do show a distinct change that would help distinguish between the canted SCI phase and its planar limit. The change 
seen in Fig. [5] clearly shows the loss of intensity for the SW modes at (0,0) after entering the planar region with 7 > 
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FIG. 6: (a) Transverse low- and high-frequency SW modes for B' = 0...2 with r\ = 7 = 4 as function of (k x /n,k y /n). The 
corresponding SF intensities over the same values are shown below. The higher frequency modes typically have a lower intensity 
in the k x direction, (b) Longitudinal modes and intensities activated by the canting of the local moments. 



8/3. Note that the transverse and longitudinal frequencies become identical in the planar phase. But for 7 < 8/3 in 
the canted phase the SW modes are cither purely transverse or longitudinal. 

Once a magnetic field is applied, the SW frequencies increase due to the enhanced stiffness of the local moment. 
Figure [6] shows the progression of the high and low SW frequencies and their corresponding SF intensities as the 
applied magnetic field increases from B' = 0...2 with 77 = 7 = 4. This transition is clearly evident in the low-frequency 
intensity at (0,0), which becomes non-zero as the field increases and the spins cant towards the z-axis. 



C. The SC2 Phase 

The angle a was given terms of the field and 7 by Eq. (|18|) . Although a does not depend on 77, changing r\ does 
modify the SW frequencies and SF intensities. The SW frequencies for this phase can be determined analytically by 
the equation given by Fishman [2l[ assuming different coefficients (Eq. (|B2|) and (|B3[) in Appendix B). As with the 
SCI, the SF intensities must be determined numerically. However, as for the FM phase, the zero- field intensities of 
the SC2 phase with a = tt/2 can be solved analytically. 

In zero field, the SW frequencies are 

w k = y^i^k ± {V + 1) (1 - cos (k x )), (24) 

where 

= (V - l)(cos(fc x ) - 1) + 2 7 (-l ± cos(^)). (25) 
The SF intensities arc given analytically by 
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(1,0) (0,0) (0,1) (1,0) (0,0) (0,1) 

(k /n,k /it) 



FIG. 7: (a) Contour plot of the low- and high-frequency SW modes for B' = with r\ = 7 = -2 as function of (k x /ir,k y /it). SF 
intensity is indicated from to highest to lowest (red - indigo). Note that the intensity of the k x is on a linear scale, while k y 
intensity is on a logarithmic scale due to a strong peak at (0,1). (b) Low- and high-frequency SW modes for B' = 0...8 with r\ 
= 7 = -2 as function of (k x / V ,k y / V) . Corresponding SF intensities over the same values. The higher-frequency modes have a 
lower intensity in the k x direction. 

Even though the angles of the local moments do not change throughout the zero-field phase, the SW frequencies and 
SF intensities depend on the interaction parameters, 7 and 77, where the zero-field intensity for each branch in Eq. 
(|2~4")) is the same. 

Figure [TJa) shows a contour plot of the SW frequencies and SF intensities for the zero-field SC2 phase with 7 = 
i] = -2. The diagram in the right-hand corner, illustrates the zero-field angles for the SC2 phase, where the local 
moments have angles a = 7r/2. The intensity of the modes in the (k x /Tr,0) region are on a linear scale, while the 
(0,fcj,/7r) intensity is set on a logarithmic scale with an intensity peak at (0,1). 

To further investigate the field-induced canting in this phase, we evaluate the SW frequencies and SF intensities 
for various magnetic field for 7 = 77 = -2 (shown in Fig. [3b)). This demonstrates the transition from the planar 
SC2 phase into a field- induced canted phase and then finally into the FM phase as a goes from ir/2 to 0. Each mode 
of the SC2 phase has both longitudinal and transverse character. This is different from the SCI phase, where the 
longitudinal modes are separate from the transverse. 

V. CONCLUSION 

In an attempt to understand the nature of competing FM and AF interactions and canted spin moments, we present 
a new technique for the determination of the static and dynamic properties for any periodic magnetic spin system. 
Using Euler angles, we determine the interactions within the local frame of reference for each spin, and apply these 
rotations to a Holstcin-Primakoff expansion to determine the classical energy and SW frequencies. A Green's function 
technique is used to determine the SF intensities for any eigcnfrcqucncy of the system. 

This technique was then applied to the generalized Villain model (GVM), which has been further generalized by 
introducing a varying interchain coupling. By studying the affects of this new interaction, we hope to gain deeper 
understanding of frustrated magnetic systems. To obtain the phase boundaries and dynamics throughout the GVM', 
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we introduce two spin configurations: SCI and SC2. 

With the phase space established, the SW frequencies and SF intensities for the different phases were determined. 
The SW frequencies were determined analytically through the whole phase space. In most cases, the SF intensities 
were determined numerically. However, we were able to provide an analytical equation for the FM and zero-field 
SC2 phase that showed how the SF intensity depends on the interactions. Using these quantities, we presented a 
cross-section of the phase space to give an overall picture of the GVM'. 

We hope that this technique will prove useful in the understanding of many different non-collincar magnetic systems. 
Magnetic heterostructures as well as frustrated AFs like CuFc02 are some of the systems that may be studied with 
this technique. 
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APPENDIX A: ROTATION COEFFICIENTS 



Using Euler rotations with angles (9 and ifj) [231 ] , each spin Si is rotated into its local frame of reference using Si 
t/.^Si, where the rotation matrix, U_ il is given by 



U; 



cos($i) cos(V>i) cos($i) sm(ipi) — sm(8i) 

- sin(V>i) eos(V>i) 

sin(#i) cos(V'i) sin(9i) sm(ipi) cos(6i) 



(Al) 



When examining the interaction between two spin operators, the overall rotation from one reference frame to another 
is 



UJT7 



pv pv 

xx xy xz 

J?ij pij J?ij 

yx yy yz 

pu pv ptj 

zx zy zz 



(A2) 



where can be obtained from Eq. (|A1[) . It should be noted that there are St different unitary matrices, one for 
each magnetic SL r. The rotation coefficients for the generalized second-order Hamiltonian (Eq. Q) are 



1 2 ^ xx ' yy \ xy r yx)) 

G r 2 s = -l(F™-F™-i(F™ + F™)). 



(A3) 



When calculating the SW intensities, one needs to multiply each element in the spin Green's function matrix by a 
rotation coefficient that describes the transverse and longitudinal contributions to the SF. The rotational coefficients 
of the transverse components (S^(O)SJ(i)) are given by 



2C 



2s-l,2r-l 



Cr 



\^xy,xy ' ^yy,yy/ > V^xx.yy ^yx,xy) ~T~ v^xy^xy ^yy,xx) 



' \ ' \ xx,yx yx,xx) \ xy,yy yy,xy) \ xx,xy ' yx,yy) \ xy,xx * yy : yx/J 



2C 



2s,2r-l 



xx,xx ' yx.yx \ xy,xy * yy,yy) \ xx.yy yx,xy) \ xy 7 xy yy.xx) 

i a'{' -r- ((~1T8 _ S~1T8 \ I Ir^iTS ' _ S^ITS i _i_ (/~1TS _j_ f^iVS \ i ir^rs _j_ r~ir s \\ 

' V ' \ xx.yx ^yXjXx) \^xy,yy ^yy,xy) ' \^xx,xy ' ^yx.yy) ' V^xy.xx ~* ^yy,yx J ) 



2C 



2s-l,2r 



(A4) 



±,=p *~^xx,xx ' ^yx,yx ' W^xy,xy ' ^yy,yy) ~T~ V^xx^yy ^yx : xy) \^xy,xy ^yy,xx) 

• ( (s~irs _ s~irs \ -j- (s-ivs ' _ r-irs \ fs-irs _|_ r-irs \ , fs^TS ■ rirs \\ 

\ ~r \^-'xx. : yx ^yx.xx) ~r \^xy,yy yy,xy) \^xx,xy ' ^yx^yy) * \^xy,xx ~* yy ,yx ) } 



2C: 



Cr 
r 



' yx t yx ' \ xy,xy ' ^yy,yy) \^x 



^yx,xy) ~r \^xy,xy ^yy,xx) 



4-?( =c (C rs — C rS \ =P (C ra — C rs \ A- (C rs -\- C rs ) — (C rs -\- c rs W 
' V i V xx,yx ^yx.xx) ~r \^xy,yy ^yy : xyJ ' v^xx.xy ' ^yx^yyJ \^xy,xx ' ^yy.yxJ J ? 
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where cd = U ab lr U cd Ls is defined in terms of the inverse of the rotation matrix elements. Note that r and s denote 
a specific SL dependent on angles 9 and tp. The rotational coefficients for the longitudinal component (S^(O)Sf (i)) 
are given by 

c\/-i1s— l,2r— 1 /~*rs r-irs \(f^ rs _i_ C TS \ 

L ^z,z ~ °13,13 °23,23 4 V°13,23 "•" °23,12,' 

n^2s,2r — 1 prs ■ /^rs i -Ir^rs i prs \ 

zo z,z — °13,13 °23,23 ' H°13,23 t °23,12/ 

(A5) 

(j/~i2s — 1.2r /irs _■_ prs i(r* TS _1_ /^ rs *\ 

zo z,z 13,13 ~ ^ 23,23 H^IS^S °23,12/' 

ZL/ z,z — °13,13 — °23, 23 t H°13,23 "+" u 23,12j- 

APPENDIX B: SPIN-WAVE FREQUENCY COEFFICIENTS 

When 7=1, the second order Hamiltonian (Eq. ((20|) ) and SW frequencies were determined analytically by Fishman 
in Ref. [21j. For the GVM' Hamiltonian with general 7, the coefficients for SCI are 

A<£' a) = 2cos(20„) +27eos(0 a -6 b ) - 2 cos 2 (6* a ) cos(k x ) + B' cos(6> a ) 

A {a,b) = A (b,a) = _ 27COS 2 (( ^ _ b)/ 2)cO8(k y ) 

A^' h) = 2 7 cos((9 a - 6 b ) - 2r)cos(26 b ) + 2?/cos 2 ((9 fc ) cos(k x ) + £'cos(0&) 
B { ^ a) =sin 2 (e a )cos(k x ) 

B (a,b) = B (b,a) = _ 27sin 2 (( ^ _ 9b) /2)cos(k y ) 

B^' b) = -ricos 2 (8 b )cos(k x ), 

while the coefficients for SC2 are given by 

^ a ' a) =2(1+7 cos(0 a - 6 b ) - cos(fc x )) + B' cos(0 a ) 

A (a, b ) = A (b,a) = _ 27COS 2 (( a _ e b) / 2 )cos(k y ) 

A$' b) = 2( 7 cos(# Q - 9 b ) - ri + V cos(k x )) + B 1 cos(9 b ) 



(Bl) 



(B2) 



B (a. b ) = B (b,a) = _ 27sin 2 (( ^ _ b)/ 2)cos(k y ) 

<< b) = 0. 



The SW frequencies are given in terms of A^'^ and B^' 8 ' by 



«3k = 4(4 a ' Q)2 + 4 M)2 - 4(4 a ' Q)2 + B^ 2 ) ) (4 a ' b)2 - 4 a ' h)2 ) + [A^ a)2 - A^ b)2 - 4(4 a ' Q)2 - B {b ' b)2 \ 



2 



+8(4 a ' a) 4 b ^ ) + 4B^ a) <' 6) ) (4 a ' b)2 + B|f' b)2 ) - 32A^' b) B^' b) (A ( ^ a) B (bM + A^' b) B^' a A 

(B4) 
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